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ABSTRACT 

The uncertainty in the geopotential model is currently one of the major error 
sources in the orbit determination of low-altitude Earth-orbiting spacecraft. 
This paper presents the results of an investigation of different geopotential 
error models and modeling approaches currently used for operational orbit 
error analysis support at the Goddard Space Flight Center (GSFC), with em- 
phasis placed on sequential orbit error analysis using a Kalman filtering algo- 
rithm. 

Several geopotential models, known as the Goddard Earth Models (GEMs), 
have been developed and used at GSFC for orbit determination. The errors m 
the geopotential models arise from the truncation errors that result from the 
omission of higher order terms (omission errors) and the errors in the spheri- 
cal harmonic coefficients themselves (commission errors). At GSFC, two error 
modeling approaches have been operationally used to analyze the effects of geo- 
potential uncertainties on the accuracy of spacecraft orbit determination - the 
lumped error modeling and uncorrelated error modeling. 

The lumped error modeling approach computes the orbit determination errors 
on the basis of either the calibrated standard deviations of a geopotential 
model’s coefficients or the weighted difference between two independently 
derived geopotential models. The uncorrelated error modeling approach treats 
the errors in the individual spherical harmonic components as uncorrelated 
error sources and computes the aggregate effect using a combination of indi- 
vidual coefficient effects. 

The study presented in this paper assesses the reasonableness of the two error 
modeling approaches in terms of global error distribution characteristics and 
orbit error analysis results. Specifically, this study presents the globa istn " 
bution of geopotential acceleration errors for several gravity error models and 
assesses the orbit determination errors resulting from these error models for 
three types of spacecraft -the Gamma Ray Observatory, the Ocean Topogra- 
phy Experiment, and the Cosmic Background Explorer. 


* This work was performed for the National Aeronautics and Space Administration (NASA)/Goddard Space 
Flight Center (GSFC) under Contract NAS 5-31500. 


6130-15 


229 



1. INTRODUCTION 


The Earth’s geopotential field, V, is commonly represented by the following spherical har- 
monic equation: 


V = — 
r 


1 + X X (Re/r) n P™ (sin 0) [Sn sin mX 

n = 2 m = 0 


+ 


Cn cos mA] 


( 1 ) 


where p 


r 


Re 


P? 


S m 

n > 


c? 


0 , X 


gravitational constant times mass of the Earth 
spacecraft orbital radius 

radius of the Earth (usually taken as the equatorial radius) 
associated Legendre function 

harmonic coefficients of degree n and order m 
geocentric latitude and east longitude 


Of computational necessity, the geopotential field represented by Equation (1) is truncated at 
some finite degree and order. Over the years, several progressively more accurate geopoten- 
tial models, known as the Goddard Earth Models (GEMs), have been developed and used at 
Goddard Space Flight Center (GSFC) for satellite orbit determination. Although the geopo- 
tential modeling accuracy has improved, it remains one of the major error sources in the orbit 
determination of low-altitude Earth-orbiting spacecraft. The errors in the geopotential force 
models arise from the truncation errors resulting from the omission of higher order terms 
(omission errors) and the errors in the spherical harmonic coefficients (commission errors). 

In deriving each GEM, a covariance matrix associated with the spherical harmonic coeffi- 
cients (or the calibrated form of this matrix) can generally be used to analyze the effects of 
commission errors on spacecraft orbit determination. However, because of the large size of 
the matrix, it is computationally expensive and, thus, is impractical for operational use. 

As an alternative, a computationally simpler modeling approach, known as the lumped error 
model (Reference 1), has been used at GSFC for orbit error analysis using the batch-weighted 
least-squares estimation method. A second approach currently being used in sequential error 
analysis, the uncorrelated error model , treats each spherical harmonic component as uncorre- 
lated and considers each to be an independent error source. 


1.1 LUMPED ERROR MODEL 


In the lumped error modeling approach (LEMA), the orbit determination errors that result 
from the geopotential model errors are computed by summing algebraically the contributions 
of the individual spherical harmonic coefficient errors. The algebraic summation assumes 
that the individual harmonic coefficient errors are fully correlated. This approach also results 
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in better computational efficiency by having to evaluate only one variational equation to ac- 
count for the geopotential errors contributed from all harmonic coefficient error sources 
(Reference 2). The error models used for this approach commonly belong to two categories: 
gravity difference and standard deviations . 

In the first category, an error model is constructed by taking the weighted difference between 
two independently derived geopotential model coefficients— the gravity difference error 
models. The rationale for such an approach is that, if two geopotential models are independ- 
ently derived, then the weighted average of the two model coefficients may be closer to the 
“truth ” and the weighted difference between the two models can be regarded as a measure o 
the error for one of the models. However, operationally, the gravity difference models have 
been usually generated based on differencing two somewhat correlated geopotential models, 
for example, GEM9-GEM7 gravity difference model. 

In the second category, the estimated errors of the geopotential coefficients -the calibrated 
standard deviations— are used as the error model. These standard deviation values are ob- 
tained from the corresponding calibrated error covariance matrix obtained for each geo- 
potential model (References 3, 4, and 5). However, it was demonstrated in References 6 and 7 
that the use of GEM9 standard deviations as a lumped error model can be faulty, because of 
the anomalous global distribution of acceleration errors that is greater in the northern hemi- 
sphere with a singularity near zero degree longitude and 60 degrees latitude. Such anomalous 
distribution is not supported, however, by other evidence (References 3 and »). 

1.2 UNCORRELATED ERROR MODEL 

The uncorrelated error modeling approach (UEM A) is generally used with the standard devi- 
ations models and assumes that each spherical harmonic coefficient can be treated as an inde- 
pendent error source. The orbital error caused by the geopotential error is then computed as 
the root sum square (RSS) of the independent error contributions from all the harmonic coef- 
ficients. This approach still requires a substantial amount of computational resources w en 
compared with the LEM A because it requires evaluating one variational equation for each o 
the spherical harmonic coefficient error sources (Reference 2). For example, for the GEMT1 
model with degree and order of 36, UEMA will require evaluating 1363 variational equa- 
tions. Although some of the coefficients are known to be correlated (Reference 8) this ap- 
proach may be reasonable (short of having to consider the correlations by including the entire 

covariance matrix). 

It was demonstrated in Reference 6, using the GEM9 standard deviations model that this 
approach produced more uniform global error distribution and more realistic orbital error 
distribution, and has, thus, been used previously in sequential error analyses performed for 

GSFC. 

The study presented in this paper investigates the LEM A and UEMA in terms of global error 
distribution characteristics and reasonableness of error magnitudes predicted for a variety o 
spacecraft orbital scenarios. The reasonableness of a geopotential error model is assessed in 
terms of global error distribution characteristics and orbit error analysis results. The geopo- 
tential error models studied are listed in Table 1 . The study results are presented in two parts. 

In the first part the global error distribution characteristics are analyzed for each of the geo- 
potentiafm^dels described in Table 1. In the second part of the study, the orbit determination 
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Table 1. Geopotential Error Models and Modeling Approaches 
Investigated in This Study 


GeopotentiaJ 

Model 


GEM 9 


Corresponding Error Models 


Model 


GEM9 Standard Deviations 
GEM9 Standard Deviations 
GEM9 - GEM7 Gravity Difference a 


Size 

(Deg. Order) 


30, 30 


Approach 


LEMA 

UEMA 

LEMA 


GEMT1 


GEMT1 Standard Deviations 
GEMT1 Standard Deviations 
GEMT1 - Clone b 


36, 36 


LEMA 

UEMA 

LEMA 


GEMT2 


GEMT2 Standard Deviations 
GEMT2 Standard Deviations 
GEMT2 - Clone b 


50, 50 


LEMA 

UEMA 

LEMA 


Notes: 


LEMA: Lumped error modeling approach 
UEMA: Uncorrelated error modeling approach 


a. GEMS - GEM7 difference model is the geopotentiai error model currently used for operational error analysis support to 

represent the 3 cvalues of GEM9 force model error. This model is derived by taking the 135% of GEM9 - GEM7 differences uo 
to degree and order of 21 to represent commission errors of the terms used in operational orbit determination and 1 00% of 
otM9 coefficients from degree and order of 22 up to 30 to represent omission errors. 


b. The clone models (References 9 and 10) are constructed such that the gravity errors computed using the difference between 
the original models (GEMT1 and GEMT2) and the respective GEMT1 and GEMT2 clone models give similar results obtained 
using the full calibrated error covariance matrices. A multiplication factor of 3 can be applied to represent 3o error values. 


errors resulting from the geopotential errors are assessed by performing linear error analysis 
using the Orbit Determination Error Analysis System (ODEAS) (Reference 2) for a variety of 
spacecraft mission scenarios. Specifically, analyses were performed for the Gamma Ray Ob- 
servatory (GRO), the Cosmic Background Explorer (COBE), and the Ocean Topography Ex- 
periment (TOPEX) types of missions, primarily in a sequential Kalman filtering mode. 

2. RESULTS 

The results of this study are presented in two parts. Section 2.1 discusses the global distribu- 
tion of geopotential acceleration errors for each of the error models listed in Table 1. Sec- 
tion 2.2 discusses the orbit determination errors resulting from geopotential errors predicted 
by different geopotential models and modeling approaches for GRO, COBE, and TOPEX 
spacecraft missions. 

2-1 GLOBAL DISTRIBUTION OF GEOPOTENTIAL ACCELERATION ERRORS 

Errors in the spherical harmonic coefficients will cause errors in computing the geopotential 
accelerations (forces) acting on an orbiting spacecraft. These acceleration errors vary as a 
function of geocentric latitude and longitude, and the error magnitude decreases with the in- 
crease in orbital height for a fixed latitude and longitude grid point in space. Figures 1 
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through 9 show the acceleration error magnitudes ( lcr values) as a function of geocentric 
latitude and east longitude at an altitude of 200 kilometers above the surface of the Earth for 
different geopotential error models and modeling approaches. The acceleration errors 
shown are in the units of mgal (10~ 5 m/sec 2 ), rounded to the nearest integer. 

For the LEMA, the acceleration errors at each of the latitude and longitude grid points are 
computed by algebraically summing the contributions from each of the sine and cosine har- 
monic coefficients at that particular location, using Equation (2). 


n max 

AaLEMA = X 
n = 2 


2 


m = 0 



AC™ 


da 

as? 



( 2 ) 


where n max represents maximum degree of the geopotential error model, and AaLEMA repre- 
sents an instantaneous geopotential acceleration error a spacecraft would experience due to 
errors in the spherical harmonic coefficients C™ and S™ , which are represented by AC n and 
AS™ in Equation (2). The acceleration error maps obtained using this modeling approach 
are shown in Figures 1 through 6. In orbit error analysis computation using the LEMA, the 
spacecraft position and velocity errors are obtained by integrating the effect of AaLEMA given 
by Equation (2) in one variational equation. 

For the UEMA, the acceleration errors are obtained by computing the RSS of the contribu- 
tions from each of the sine and cosine harmonic coefficients rather than summing them 
algebraically. This approach is shown mathematically in Equation (3) where AayEMA repre- 
sents the acceleration error magnitude obtained using the UEMA. The acceleration error 
maps obtained using this modeling approach are shown in Figures 7 through 9. 


AauEMA 



da 

dO 1 


AC„ 


da 

dSS* 



( 3 ) 


It should be noted that in orbit error analysis computation using the UEMA, the spacecraft 
position and velocity errors are obtained by integrating the individual vectors in the summa- 
tion (i.e., — AC? and AS{? ) in separate variational equations and computing the 
dC™ dS n 

RSS of all individual error contributions. 

It is interesting to see from Figures 1 through 3 that all acceleration error maps obtained for 
the LEMA using the GEM9, GEMT1, and GEMT2 standard deviations models tend to show 
the same anomalous error distribution (i.e., a high concentration of large errors localized in a 
particular region around zero degree longitude in the northern hemisphere and the rest of the 
globe relatively error-free). The world map of gravity anomalies computed from GEM9 and 
GEM10 models (Reference 3), the gravity anomalies computed from GEM7 and GEM8, and 
the gravimetry data (Reference 8) do not support such anomalous distribution. This implies 
that the orbit determination errors predicted by using such error models will also be anoma- 
lous. 
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Figure 4. Gravitational Acceleration Errors (mgal) Based on the LEMA Using 
the GEM9-GEM7 Difference Error Model 
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the GEMT1 -Clone Difference Error Model 
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Figure 7. Gravitational Acceleration Errors (mgai) Based 
the GEM9 Standard Deviations Model 


280 300 320 340 360 

on the UEMA Using 


<LAT> 


90 

80 

70 

60 

50 

40 

30 

20 

10 

0 

-10 

•20 

-30 

-40 

•50 

•60 

•70 

•SO 

-90 


3 3 
3 3 


22222222 


3 

4 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4 

3 

2 


2 2 


J-— 1-— i— 1— 




- I— i l I- 


, I— I— I_ 


<LOXG> 


60 


80 


100 120 160 ISO 200 220 240 260 280 300 320 340 360 


Figure 8. 


<LAT> 


Gravitational Acceleration Errors (mgai) Based on the UEMA Using 
the GEMT1 Standard Deviations Model 


90 

SO 

70 

60 

50 

40 

30 

20 

10 

0 

• 10 
-20 
•30 
-40 
•50 
-60 
-70 
-80 
-90 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


3 3 
3 3 


3 3 
3 3 


3 3 
3 3 


3 


3 
3 
3 
3 

3 3 

3 3 

3 3 

3 3 

3 3 

3 3 3 

3 3 3 

3 3 

3 3 

3 3 

3 3 


3 } 


3 

3 

4 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 

4 

3 

3 


3 

3 

4 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 

4 

3 

3 


3 

3 

4 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 

4 

3 

3 


3 

3 

4 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 3 

3 3 

3 3 

4 3 

3 3 

3 3 


' I— i — i — I — j — ] — i — ; — j — , — , — j — j — j — | j.. 


■ I— - I— - I-- I-— !- 


<LONC» 

C 


40 60 


300 . J8C ZX 320 2*0 260 


Figure 9. 


0 300 320 340 360 

Gravitational Acceleration Errors (mgai) Based on the UEMA Using 
the GEMT2 Standard Deviations Model 


6130-15 


236 



On the other hand, the acceleration error distributions obtained for the UEMA seem to be 
too uniform and symmetric to be realistic. (See Figures 7 through 9.) Comparisons between 
GEM models and surface gravity measurement data (References 3 and 8) show that a more 
realistic error distribution may be somewhere between these two cases. 

It is also interesting to see that the LEMA based on the gravity difference models (i.e., 
GEM9-GEM7, GEMTl-Clone, and GEMT2-Clone models) give rise to error maps showing 
random distributions in error magnitudes, (See Figures 4 through 6.) In the absence of the 
absolute standard, the gravity difference models appear to be more realistic in terms of global 
acceleration errors distribution. 

2.2 ORBIT DETERMINATION ERRORS RESULTING FROM GEOPOTENTIAL 
MODELING ERRORS 

To assess the orbit determination errors that result from the geopotential modeling errors, 
linear error analyses were performed using ODEAS for a variety of spacecraft missions with 
different altitudes and inclinations. Specifically, error analyses were performed for GRO, 
COBE, and TOPEX missions assuming an orbit determination mode using a sequential 
Kalman filter. Batch-weighted least-squares results are also shown for COBE spacecraft. 

The orbital elements at epoch for GRO, COBE, and TOPEX, their respective spacecraft pa- 
rameters, force models, and the integration parameters used in the study are shown in Table 2. 
Tracking simulations were performed for each spacecraft assuming a Tracking and Data Relay 
Satellite System (TDRSS) tracking scenario, whereby the user spacecraft (GRO, COBE, or 
TOPEX) is tracked for 2 days via TDRS-East (41 degrees west longitude) and TDRS-West 
(171 degrees west longitude) with tracking pass lengths of 20 minutes per user spacecraft rev- 
olution. Measurements include one range and one range rate every 10 seconds. 

Error analyses were performed by estimating the orbital elements of the user spacecraft and 
the atmospheric drag scaling parameter. Because the primary interest of this study is in 
orbital errors resulting from the geopotential error sources, only the characteristics of geopo- 
tential errors on orbit determination accuracy are presented. The orbital errors presented 
below are based on 3a errors. 

Figures 10a through 10c illustrate the GRO (altitude: 450 km) spacecraft orbit determina- 
tion errors resulting from the GEM9, GEMT1, and GEMT2 geopotential model errors, re- 
spectively. Figures 11 and 12 illustrate similar plots of orbit determination errors for COBE 
(altitude: 900 km) and TOPEX (altitude: 1340 km) spacecraft, respectively. Table 3 shows 
the root mean square (RMS) values, standard deviations, and maximum values of the orbital 
errors after the initial transient period (taken to be 6 hours past epoch). 

The RMS and maximum errors resulting from a given geopotential error model are expected 
to decrease with increasing orbital altitudes from GRO to TOPEX. The orbital errors for a 
given spacecraft are also expected to decrease from GEM9 to the more refined 
models— GEMT1 and GEMT2. The observations made for each of the geopotential model- 
ing approaches are summarized below. 

1. UEMA results using the standard deviations models generally produce rather uni- 
form error distributions at the expense of substantially higher computer processing times 
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Table 2. Spacecraft Orbital Elements and Parameters 



GRO 

COBE 

TOPEX 

EPOCH: (YYMMDD) 

920201 

921109 

920901 

(HHMMSS) 

000000 

000000 

000000 

ORBITAL ELEMENTS AT EPOCH 

Semi-major axis (km) 

6828.15 

7278.047 

7718.756 

Eccentricity 

0.0001 

0.0003238 

0.00112783 

Inclination (deg) 

28.5 

99.0296 

63.14283 

R. A. of ascending node (deg) 

0.0 

86.8450 

0.0 

Argument of perigee (deg) 

0.0 

94.9595 

0.0 

Mean anomaly (deg) 

0.0 

270.8221 

0.0 

Orbital Period (minutes) 

93.59 

102.99 

112.48 

SPACECRAFT PARAMETERS 

Drag coefficient (Cq) 

2.2 

2.3 

2.3 

Spacecraft area to mass 

0.003347 

0.006487 

0.006487 

(m 2 /kg) 

Solar reflectivity (C R ) 

1.5 

1.3 

1.3 

FORCE MODELING 

GEM9 30 X 30 for reference trajectory computation 
Solar radiation pressure 
Solid Earth tide 

Atmospheric drag force using Harris-Priester density model 
(Mean F 10 7 solar flux level - 1 75 X 1 0 -22 watts/m 2 /Hz) 




ORBIT INTEGRATION 


Integrator: Fourth-order Runge-Kutta Integrator 
Step size: 60 seconds 


TRACKING SCENARIO 

20 min/user spacecraft orbit from TDRS-East or TDRS-West 

Range noise - 3.0 m 

Range-rate noise - 0.00282 m/sec 

Date rate - 1 range and range rate per 10 seconds 

For batch-weighted least-squares method, the data weights are set to be the same as for data noise 
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Figure 10b. GEMT1 Geopotential Errors for GRO Processed in Sequential Mode 
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Figure 11a. GEM9 Geopotential Errors for COBE Processed in Sequential Mode 
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Figure 11c. GEMT2 Geopotential Errors for COBE Processed in Sequential Mode 
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Figure 12a. GEM9 Geopotential Errors for TOPEX Processed in Sequential Mode 



--- GEMT1 SIGMAS LEMA GEMTl SIGMAS UEMA — - GEMT 1 -CLONE LEMA | 


Figure 12b. GEMTl Geopotential Errors for TOPEX Processed in Sequential Mode 
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Figure 12c. GEMT2 Geopotential Errors for TOPEX Processed in Sequential Mode 
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Table 3. Steady-State RMSs, Maximums, and Standard Deviations 
of Orbital Errors in Sequential Filtering Mode 



LEMAL 

C 

JSING STANDARD 
EVIATIONS 

LEMA USING GRAVITY 
DIFFERENCES 

UEMA USING STANDARD 
DEVIATIONS 


GEM9 

GEMT1 

GEMT2 

GEM9 a 

GEMT1 b 

GEMT2 C 

GEM9 

GEMT1 

GEMT2 

RMS 










GRO 

20.462 

8.578 

5.693 

85.155 

105.188 

96,495 

119.717 

83.700 

76.112 

COBE 

96.004 

30.581 

23.799 

56.067 

26.588 

16.786 

48.529 

25.288 

19.211 

TOPEX 

6.59 a 

3.463 

3.027 

22.459 

3.527 

1.813 

22.111 

9.500 

7.689 

Maximum 










GRO 

42.344 

23.146 

15,784 

200.430 

288.670 

235.400 

227.350 

135.420 

117.960 

COBE 

323.620 

86.041 

66.670 

181.480 

60.767 

53.962 

85.404 

43.749 

34.032 

TOPEX 

22.229 

8.848 

9.052 

63.587 

10.399 

6.866 

42.564 

18.198 

15.028 

Standard 










Deviation 










GRO 

8.462 

4.894 

3.148 

33.998 

56.350 

40.462 

22.190 

14.929 

14.059 ; 

COBE 

50.265 

16.092 

12.824 

26.493 

14.655 

9.964 

15.127 

7.129 

5.356 < 

TOPEX 

3.498 

1.841 

1.702 

12.737 

1.848 

1.199 

5.498 

2.274 

1.881 < 


Units are in meters 


Nate s; 


a. GEM9-GEM7 difference model is used to represent GEM9 errors 

b. GEMT1 -Clone difference model Is used to represent GEMT1 errors 

c. GEMT2-Clone difference model is used to represent GEMT2 errors 


(generally requiring three to seven times as much central processing unit (CPU) time as the 
corresponding LEMA results). The uniformity is characterized by the relatively small stand- 
ard deviation values when compared with the corresponding RMS values for all the spacecraft 
studied. Tkble 3 also indicates that the orbital errors predicted by the UEMA follow the ex- 
pected trends. 

2. LEMA results using standard deviations models produce either large error spikes (as 
in the case of COBE) or relatively small errors (as in the case of GRO) when compared with 
the corresponding UEMA results. These results are expected because of the anomalous dis- 
tributions of geopotential acceleration errors as illustrated in Figures 1 through 3. The effect 
of such distribution is that the orbital errors predicted for spacecraft with high orbital inclina- 
tion, such as COBE (99 degrees), will result in large local error spikes. In contrast, the orbital 
errors predicted for spacecraft with lower orbital inclinations, such as GRO (28.5 degrees), 
can be unrealistically small. Such error characteristics are not supported by operational expe- 
rience. 

3. LEMA results using the gravity difference models are found to exhibit different be- 
haviors for different spacecraft orbital scenarios. From Table 3, it can be determined that the 
RMS and the maximum errors resulting from a given geopotential error model decrease, as 
expected, when going from a low-altitude orbit of GRO to a high-altitude orbit of TOPEX. 
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However, the orbital errors do not follow the expected decreasing trend from GEM9 to the 
more refined GEMT1 and GEMT2 error models for GRO, where the orbital errors resulting 
from the GEMTl-Clone and GEMT2-Clone error models are larger than those of 
GEM9-GEM7 error model. The error analysis results specific to each of the gravity differ- 
ence models are summarized below: 

• GEM9-GEM7 difference model produces reasonably close agreement with the 
GEM9 UEMA results for GRO and COBE spacecraft. For TOPEX spacecraft, 
while GEM9-GEM7 produces RMS orbital errors close to the UEMA results 
(Table 3), relatively large and spiky errors are observed in orbital error distribu- 
tions. (See Figures 10a, 11a, and 12a.) 

• GEMTl-Clone model produces results comparable with the GEMT1 UEMA re- 
sults for COBE but produces comparatively higher errors in GRO and compara- 
tively lower errors in TOPEX. (See Figures 10b, lib, and 12b.) 

• GEMT2-Clone model produces results comparable with the GEMT2 UEMA re- 
sults for COBE but also produces comparatively higher errors in GRO and lower 
errors in TOPEX. (See Figures 10c, 11c, and 12c.) 

To assess the reasonableness of the error magnitudes predicted by each error modeling ap- 
proach, the orbital errors must be compared with independent results. Reference 5 provides 
independent error analysis results for TOPEX-type spacecraft for GEMT1 and GEMT2 
models. These results were obtained using the first-order analytical perturbation theory of 
Kaula (Reference 1 1) and fully calibrated covariance matrices of GEMT1 and GEMT2 solu- 
tions, assuming a 10-day TOPEX arc. 

Because such analyses do not account for the effect of imperfect tracking scenarios, the 
orbital errors thus obtained can be regarded as errors due entirely to the geopotential model 
under a perfect tracking condition. According to Reference 5, the projected orbital errors 
resulting from GEMT1 and GEMT2 model errors for a 10-day TOPEX arc are 1.9 and 
0.9 meters (lo) , respectively. This translates to 5.7 and 2.7 meters (3a) as shown in Table 4. 

For comparison purposes, additional TOPEX error analysis was performed using the 
ODEAS sequential processing mode assuming a continuous tracking of TOPEX from 
TDRS-East and TDRS-West for 2 days using range and range-rate data. Continuous tracking 
was used to minimize the effect of imperfect tracking conditions on orbital errors. Table 4 
lists the RMS orbital errors obtained after the initial transient period using different GEMT1 
and GEMT2 error models. 

The results indicate that, for TOPEX-type spacecraft, the projected orbital errors quoted in 
Reference 5 lie in between the UEMA and the LEMA results with UEMA results being on 
the high side. 

Figures 13a through 13c show the batch-weighted least-squares orbit determination errors 
for COBE spacecraft resulting from GEM9, GEMT1, and GEMT2 error models, respec- 
tively. Error analysis results are presented assuming a 2-day definitive data for all cases. 
Again, it can be observed that the LEMA using standard deviations models produces very 
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Table 4. Comparison of TOPEX Orbital Errors (3a values) Using Different 
Geopotential Error Modeling Approaches 


Geopotential 

Model 

Projected orbital 
errors for 1 0-day 
arc (meters) 
(Reference 5) 

Steady-state RMS 
orbital errors for a 2-day 
continuous tracking (meters) 

UEMA 

LEMA 3 

GEMT1 

5.7 

7.4 

3.2 

GEMT2 

2.7 

6.1 

1.3 


Motes; 

a. GEMT1 -Clone and GEMT2-Cione gravity difference models are used to 
represent geopotential errors due to GEMT1 and GEMT2 models, 
respectively. 

large spikes in orbit errors. However, the LEMA using the gravity difference models 
produces results similar to the corresponding UEMA results. This was also found to be the 
case for COBE spacecraft in sequential processing mode. (See Figures 11a through 11c.) 

Figures 14a through 14c compare COBE orbit determination errors processed in sequential 
Kalman filtering with batch-weighted least-squares modes using identical tracking schedules 
in both modes. As expected, the effects of local error spikes are less severe in the 
batch-weighted least-squares mode than in the sequential Kalman filtering mode, and the 
batch solutions agree with the corresponding sequential solutions at the end of the data arc. 

3. CONCLUSIONS AND RECOMMENDATIONS 

From the observation of geopotential acceleration error maps, two primary conclusions can 
be drawn. 

1. The LEMA using standard deviations models produces an anomalous trend having 
a high concentration of large acceleration errors that are localized in a particular region 
around zero degree longitude in the northern hemisphere and the rest of the globe remains 
relatively error-free. Such an error distribution can give rise to spurious local error spikes for 
high-inclination spacecraft (such as COBE) and unrealistically small errors for low- 
inclination spacecraft (such as GRO). It can, therefore, be concluded that the standard devi- 
ations models should not be used in conjunction with the LEMA. 

2. However, the UEMA using standard deviations models is found to produce uni- 
form acceleration error bands symmetric in northern and southern hemispheres. This may, 
again, be unrealistic. On the other hand, the LEMA using gravity difference models gives 
error maps showing random distributions in error magnitudes. In the absence of the absolute 
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Figure 13a. GEM9 Geopotential Errors for COBE Processed in Batch Mode 
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Figure 14b. Geopotential Errors for COBE Using GEMT1 -Clone Model 
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Figure 14c. Geopotential Errors for COBE Using GEMT2-Clone Model 
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standard, the error characteristics shown in the gravity difference models appear to be the 
more reasonable representation of geopotential errors. 

From the error analysis results, three additional conclusions can be drawn. 

1. The geopotential error models that are based on the UEMA generally produce 
fairly uniform error distributions usually free from unrealistic error spikes at the expense of 
substantially higher computer processing times (generally requiring three to seven times as 
much CPU as the LEMA). The uniformity is characterized by relatively small standard devi- 
ation values when compared with the corresponding RMS values. This was found to be the 
case in both sequential and batch processing modes. Also, the orbital errors predicted by this 
modeling approach generally follow the expected trends— meaning that, for a given geopo- 
tential error model, the RMS and the maximum orbital errors decrease with increasing 
orbital altitudes from GRO to TOPEX; and for a given spacecraft, the orbital errors decrease 
from GEM9 to the more refined models (GEMT1 and GEMT2). 

2. The difference error models in LEMA generally produce higher fluctuation in er- 
rors than the UEMA. This is characterized by relatively large standard deviation values when 
compared with the corresponding RMS values. 

3. Generally, the error characteristics produced by different geopotential error mod- 
eling approaches are not sensitive to orbit error estimation mode. For COBE spacecraft, the 
error signatures produced by the batch-weighted least-squares mode are similar to those of 
sequential Kalman filtering mode, except that the error magnitudes and the error fluctuations 
are smaller in batch mode. This is expected because of the data-smoothing effect realized in 
the batch estimation mode. 

It is difficult to suggest a general error modeling approach that can be used for all spacecraft 
scenarios. As discussed earlier, if central processing unit (CPU) and computing resources are 
not of a concern, error analysis should be performed using the entire covariance matrix of 
spherical harmonic coefficients, to account correctly for the correlations that exist among 
these coefficients. In the absence of this, the UEMA may be a reasonable error modeling 
approach for most spacecraft scenarios because the orbital errors predicted by UEMA are 
mostly free from unrealistic error spikes and follow the expected trends as discussed earlier. 
The disadvantage of this approach is still the substantially higher CPU requirement over the 
LEMA for operational error analysis support. 

The LEMA using difference models, while computationally efficient, is found to exhibit dif- 
ferent characteristics for different spacecraft. GEM9-GEM7 seems to be a good error model 
to represent GEM9 modeling errors for GRO and COBE spacecraft, but it is not very reason- 
able for TOPEX because of some sporadic, large error spikes. GEMTl-Clone and 
GEMT2-Clone models seem to be good error models for representing GEMT1 and GEMT2 
errors, respectively, for TOPEX and COBE spacecraft, but tend to produce large error spikes 
in GRO spacecraft. 

As discussed previously, geopotential error models based on either the LEMA or the UEMA 
do not properly take into account the correlations that exist among the harmonic coefficients. 
If computer resources permit, future studies should include orbital error analyses using the 
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entire covariance matrix of spherical harmonic coefficients to take proper account of the cor- 
relations. This will also provide a reference point for calibrating other geopotential error 
modeling approaches. Future study should also include the development of an improved 
error modeling approach feasible for operational use. Such a modeling approach will be 
based on a certain combination of LEMA and UEMA with proper correlation information 
obtained from the calibrated error covariance. To augment the error analysis results, the au- 
thors also recommend performing sequential and batch orbit determination studies using 
real tracking data. 
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